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Abstract 

We study non-linear primordial adiabatic and isocurvature perturbations and their non-Gaussianity. 
After giving a general formulation in the context of an extended 5N formalism, we analyse in detail 
two illustrative examples. The first is a mixed curvaton-inflaton scenario in which fluctuations of both 
the inflaton and a curvaton (a light isocurvature field during inflation) contribute to the primordial 
density perturbation. The second example is that of double inflation involving two decoupled massive 
scalar fields during inflation. In the mixed curvaton-inflaton scenario we find that the bispectrum of 
primordial isocurvature perturbations may be large and comparable to the bispectrum of adiabatic 
curvature perturbations. 



1 Introduction 

With recent cosmic microwave background (CMB) anisotropy data due to the WMAP satellite 
and the further improved data expected from the Planck satellite, our knowledge of primordial 
cosmological perturbations is becoming more and more precise. This influx of data has stimulated 
the study of models whose predictions differ from the simplest models of single field slow-roll inflation. 
To discriminate between these models, a particularly important observable is the amplitude (and the 
shape) of non-Gaussianity of the CMB anisotropies. Another crucial property, potentially observable 
in the CMB data, would be the presence of a primordial isocurvature (or entropy) component as it 
would require a multi- field scenario for the origin of the primordial fluctuations. 

The purpose of this paper is to investigate the impact of non-adiabatic fluctuations during in- 
flation on the predicted non-Gaussianity of primordial density perturbations, including primordial 



isocurvature matter perturbations as well as adiabatic modes which would contribute to the bis- 
pectrum and higher-order correlations in the CMB anisotropics. Isocurvature perturbations could 
have large departures from Gaussianity while remaining sub-dominant in the linear perturbation 
spectrum [1, 2]. 

To study primordial non-Gaussianity, one needs to study relativistic cosmological perturbations 
beyond linear order and there has been considerable progress in this field in recent years. On scales 
larger than the Hubble radius, the non-linear evolution of perturbations generated during inflation is 
compactly described in terms of the perturbed expansion from an initial hypersurface (usually taken 
at Hubble crossing during inflation) up to a final uniform-density hypersurface (usually during the 
radiation-dominated era) - the so-called (5A^-formalism [3]. This is particularly useful for evaluating 
the primordial non-Gaussianity generated on large scales [4]. 

As we show in this paper, one can easily extend the (5iV-formalism to describe the non-Gaussiani- 
ties of the non-linearly evolved primordial perturbations including isocurvature fluctuations. In order 
to illustrate our general, but formal, result, we study two emblematic examples of multi-field sce- 
narios, which can generate isocurvature fluctuations in addition to the usual adiabatic fluctuations. 
The first example is the curvaton scenario [1, 5]. Previous works, e.g. [6, 7, 8], have investigated 
non-Gaussianity and isocurvature perturbations in this scenario, but in our case, we do not assume 
that the contribution of inflaton fluctuations to the CMB anisotropics is negligible. In this so-called 
mixed inflaton-curvaton setup [9, 10], the isocurvature mode is not necessarily constrained by the 
data to be zero, in contrast with the conclusion of [8]. Our second example is a model of double 
inflation [11] with two uncoupled massive scalar fields that drive in turn inflation. In contrast with 
the previous example, the final isocurvature perturbation depends on both scalar field fluctuations 
during inflation, but it can still be determined analytically at second order. 

The adiabatic and isocurvature perturbations we refer to above correspond to the primordial adi- 
abatic and isocurvature perturbations defined during the standard radiation era, i.e., after inflation 
and after the curvaton decay, if any. These perturbations can be related, but are not equivalent, to 
the instantaneous adiabatic and isocurvature (or entropy) field perturbations which can be defined 
during inflation by decomposing the perturbations along the directions, respectively, parallel and 
orthogonal to the inflationary trajectory in field space (see [12, 13] in the linear case and [14, 15] in 
the non-linear case) . The instantaneous isocurvature perturbation during inflation is not necessarily 
converted into an isocurvature perturbation after inflation. However, even if the isocurvature per- 
turbation during inflation does not survive, it can have a strong impact on the resulting primordial 
adiabatic perturbation and its non-Gaussianity, as illustrated, for instance, recently in the context 
of multi-field Dirac-Born-Infeld inflation [16]. 

The outline of the paper is the following. In the next section, we introduce the non-linear 
definitions of the primordial adiabatic and isocurvature perturbations and show how they are related 
to the primordial scalar field fluctuations in a very general multi-field inflation framework. The 
following section is devoted to the study of the mixed eurvaton-inflaton scenario. We then consider, 
in the fourth section, the ease of double inflation with two decoupled massive scalar fields. We discuss 
our results in the final section. In Appendix A we give some details of the calculations of section 
IV, while in Appendix B we review the decomposition into the adiabatic and entropy components of 
the field perturbations and their equations of motion at second order, and we compute their 3-point 
correlation functions. Finally, in the last appendix we give general expressions obtained using the 
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(JAT-formalism for the primordial power spectra and bispectra at leading order. 

While this paper was being written up, similar results have been obtained by Kawasaki et al. who 
use the JiV-approach to calculate primordial non-Gaussianity of isocurvature perturbations and in 
particular axion isocurvature perturbations [17] and baryon isocurvature perturbations [18]. 



2 Adiabatic and isocurvature perturbations 

A powerful technique to compute the non-linear primordial perturbations on large scales is the 
(JAT-formalism [3, 4]. The idea is to use solutions to the homogeneous FRW cosmology in order to 
calculate the integrated expansion on large scales from some initial state to a final state of fixed 
energy density. 

The 5A^-formalism is closely related to the notion of non-linear curvature perturbation on uniform 
density hypersurfaces, which can be defined in a geometrical and covariant way as shown in [19, 20]. 
Indeed, in the case of a perfect fluid characterized by the energy density p, the pressure P and the 
four- velocity u", the conservation law for the energy-momentum tensor. 



implies that the covector 



satisfies the relation 



where we have defined 



Van = 0, Tab = {P + P) UaUb + PQab, (1) 



C„ = Vaa - jVaP (2) 
a-r„a = -^(VaP-|Vap), (3) 



e = Vau", a 



IJdre, (4) 



and where a dot denotes a Lie derivative along n", which is equivalent to an ordinary derivative for 
scalar quantities (e.g. p = u"'V aP)- This result is valid for any spacetime geometry and does not 
depend on Einstein's equations. In the cosmological context, a can be interpreted as a non-linear 
generalization, according to an observer following the fiuid, of the number of e-folds of the scale 
factor. 

The covector C,a can be defined for the global cosmological fluid or for any of the individual cos- 
mological fluids (the case of interacting fluids is discussed in [21]). Using the non-linear conservation 
equation 

/j = -3d(p + P), (5) 
which follows from u'''VaT\ = 0, one can re-express Q in the form 

If w = P/p is constant, the above covector is a total gradient and can be written as 



Ca = V„ 



1 

3{l + w) 



(7) 
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On scales larger than the Hubble radius, the above definitions are equivalent to the non-linear 
curvature perturbation on uniform density hypersurfaces as defined in [22, 23], 

where N = a and H = a = d/a is the Hubble rate of the Friedmann metric ds'^ = —dt^ + a^{t)dx^ . 
The above equation is simply the integrated version of (2), or of (6). 

In the following, we will be mainly interested in non-linear isocurvature, or entropy, perturba- 
tions. For simplicity, we will consider only cold dark matter (CDM) isocurvature perturbations and 
assume that the Universe, in the standard eras, is filled with only two fluids: the radiation fluid and 
the CDM fluid. Our analysis can be easily extended to other types of isocurvature perturbations. 

It will be useful to distinguish the non-linear curvature perturbation C, of the total fluid, which 
describes the primordial adiabatic perturbation, from the non-linear perturbations Q and Cm de- 
scribing respectively the radiation fluid {wj. = 1/3) and the cold dark matter (CDM) fluid {wj^ = 0), 
which are given, according to our definitions (7) or (8), by 

C. = .iV + iln(|), (9) 
U = SN + l\n(^) , (10) 



3 VPri 

where a bar denotes a homogeneous quantity. 

In the radiation dominated era, the adiabatic perturbation coincides with ^r, whereas the CDM 
isocurvature perturbation is characterized by the non-linear perturbation 

which can be expanded in terms of the density contrasts Sj- = Sp^/pT and S^a = Sp^/pjn, 

S^ = 6^-ls.-^6l + l6',+... (12) 

Note that these expressions are independent of the hypersurface on which the density perturbations 
are defined. 

Since our goal is to relate the perturbations in the radiation era to the perturbations produced 
during an inflationary era, it is important to generalize Eq. (8) for scalar flclds. In this case a 
convenient description is in terms of the (relative) comoving curvature perturbation,^ 

nA = SN- Th^, (13) 

which is the curvature perturbation on constant ipA hypersurface. In slow-roll inflation, the initial 

state of the system - when the cosmological perturbations are produced - is defined only by the 
scalar field values, if a*, on an initial spatially-flat hypersurface, where with a star we denote that 



^Note that the convention adopted here is that TZa has the same sign as (a, such that, in the single field 
case, TZ = ( on large scales. 
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we evaluate the quantity at Hubble crossing k = aH. One can then calculate the number of e-folds, 
or integrated expansion, N^'^^\ from this initial state to a "final" hypersurface characterized by 
the "final" scalar field amplitudes (pA- By choosing the final hypersurface to be of imiform A-field, 
one can write TZa as a perturbative expansion in terms of the initial field fluctuations 5(pA*, whose 
correlation properties must be known. Equation (13) thus becomes 

Ua = <5iV('^^) = N^X^hifA* + ^N^Xb^'PaJ'Pb* + ... , (14) 

where N^^^^ = dN^'^^^ /d(pA*, etc. This is a particular application of the (5iV-formalism that gen- 
eralizes the usual expansion of N defined on a final total uniform density hypersurface. Note that 
when there are several scalar fields, TZa can be different from the relative curvature perturbation on 
uniform density hypersurfaces (a- Indeed, the uniform density and uniform field hypersurfaces do 
not always coincide even on large scales [24]. 

However, the total comoving and uniform density hypersurfaces coincide on large scales at second 
[25] and non-linear order [22, 15] and ( is generally used to describe the adiabatic perturbation also 
for scalar fields. The curvature perturbation on uniform density hypersurfaces will be given now 
as the standard perturbative expansion of N defined on a final uniform density hypersurface. Thus 
one can rewrite ^ in terms of the expansion [4] 

( = SN = N^aS(Pa* + ^N^abS(Pa*S(Pb* + ... ■ (15) 

Similarly, the non- linear isocurvature perturbation (11) can be given in terms of the difference in 
the non-linear expansion, = 36{AN), where AN = N^^^ — N^^\ between final hypersurfaces of 
uniform matter density and uniform radiation density 

5m = 3 ((5iV(°^) - 6N^'^^ = SANaS^a* + ^ANabS^^a^S^b* + • • • ■ (16) 

In the following sections we apply these definitions to two examples: the curvaton model and 
double inflation. Although the previous expressions hold non-linearly, we will concentrate on a 
second order expansion, which is expected to give the leading order terms for the 3-point correla- 
tion properties. We will assume that the initial field perturbations (on scales close to the horizon 
scale during inflation) are independent, Gaussian random flelds. Thus any non-Gaussianity of the 
curvature perturbations will arise from the non-linear terms in Eqs. (15) and (16). This is a good ap- 
proximation for weakly-coupled scalar flelds (with canonical kinetic terms) during slow-roll inflation 
[26] but may break-down for scalar flelds with non-standard kinetic terms. 

3 Mixed inflaton and curvaton perturbations 

As a flrst application of the general formalism presented in the previous section, we consider a 
curvaton scenario [5], or more precisely a mixed inflaton and curvaton scenario [9, 10] as we will take 
into account both the perturbations generated by the inflaton fleld driving inflation and the curvaton. 
The curvaton is a weakly coupled scalar field, %, which is light relative to the Hubble rate during 
infiation, and hence acquires an almost scale-invariant spectrum and effectively Gaussian distribution 
of perturbations, Sx, during infiation. After inflation the Hubble rate drops and eventually the 
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curvaton becomes non-relativistic so that its energy density grows relative to radiation, until it 
contributes a significant fraction of the total energy density, Q,^ = Px/p, before it decays. Hence the 
initial curvaton field perturbations on large scales can give rise to a primordial density perturbation 
after it decays. 

The non-relativistic curvaton (mass m 3> H), before it decays, can be described by a pressureless, 
non-interacting fiuid with energy density 

Px = ' (17) 

where % is the rms amplitude of the curvaton field, which oscillates on a timescale m"^ much less 
than the Hubble time . Making use of Eq. (10) for the oscillating curvaton to rewrite its local 
density in terms of its homogeneous value and the inhomogeneous expansion perturbation, ^AT, we 
have 

= p,e3(Cx-^A^) . (18) 

In the post-inflation era where the curvaton is still subdominant, the spatially flat hypersurfaces are 
characterized by bN = (^inf, where Cinf corresponds to the adiabatic perturbation generated by the 
inflaton fiuctuations. On such a hypersurface, the curvaton energy density can be written as 

p^e3(Cx-Ci„f) = p^^S^ = (x + 5xf ■ (19) 

where Sy. = 3((^^ — C,\^i) is the entropy perturbation of the curvaton. 

As long as the curvaton is subdominant and weakly-coupled, so that we may neglect self- 
interactions, the evolution equation for x is linear, and the field perturbation 5x obeys the same 
evolution equation on super-Hubble scales as the background expectation value x- this case it 
is well known that the ratio 5x/x remains unchanged as long as the curvaton is subdominant [11]. 
This result holds also at second order in the perturbation 5x, as shown in [15] and in Appendix B, 
where we have written the evolution equation of an entropy field perturbation. Thus, expanding 
Eq. (19) at second order we obtain 



Note that we will assume that the initial curvaton field perturbations, 6x*, are strictly Gaussian, as 
would be expected for a weakly coupled field. 

The precise density perturbation produced after the curvaton decays can be calculated numeri- 
cally [27, 28, 23], but it can also be estimated analytically using the sudden-decay approximation [6], 
which assumes that the curvaton decays suddenly on a spatial hypersurface of uniform total energy 
density. Any initial inflaton perturbation gives rise to a perturbation in the radiation energy density 
before the curvaton decay, which we denote by pR. Similarly to Eqs. (9) and (10) we can write 

PR = p-«e^(^--^^) , Px = p-,e^(^-^^) . (21) 

On the decay hypersurface characterized by Px + PR = Pr and thus 5N = where is the total 
curvature perturbation after the decay, we find [23] 

^X)decay6 '''^^ ~ ^X)<iecay)6 ^^'"^ ^''•^ = 1 , (22) 
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where f^^.decay = PxKPx + Pr)- 

Expanding Eq. (22) at first order we obtain 



Cr = rCx + (l-r)Cinf, (23) 
where r = 3nj^_decay/(4 — J^)x,decay Up to second-order we obtain 

C = ^C. + (1 - r)C,., + '•'^-'•f -^'•^ - = On, + + ^<^-;f + '-' 5; . (24) 

The entropy perturbation (20) contains a linear part Sq which is Gaussian and a second order part 
which is quadratic in Sq- 

S^ = Sg- \sl , where Sq = 2^ . (25) 

Substituting in (24) we then have 

Cr = Cinf + ^SG + ^Q-2r-r2^5S. (26) 

Keeping only the linear part of the above relation, one finds that the power spectrum for the 
primordial adiabatic perturbation (^r can be expressed as 



where the entropy power spectrum amplitude is given by 



'PSa = -^[l6] ■ (28) 



1 fHA^ 



In the case of single field inflation we have 

where e^, = —H^/H^ is the usual slow-roll parameter during inflation and Mp = (SvrG)^^ is the 
reduced Planck mass. In order to compare the relative contributions of the inflaton and of the 
curvaton in the final power spectrum (27), it is useful to introduce the dimensionless parameter [29] 

so that P^j. = (1 + A)'P^.^j. If A ^ 1, one recovers the standard curvaton scenario where the inflaton 
perturbations can be ignored: since r and are bounded by 1, this requires x* ^ Mp. A value of 
A of order 1 or smaller is possible if r or e* are sufficiently small and/or x* is of the order of Mp. 
In the present work, we will always assume x* ^ Mp. If this is not the case the curvaton starts 
to oscillate at about the same time as it decays and cannot be described as a dust field (see [9] for 
details) . 
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In slow-roll inflation the 3-point function of the inflaton perturbations, Cmf , is suppressed by slow- 
roll parameters [30, 31] and large non-Gaussianities can arise only from the curvaton contribution. 
Indeed, the 3-point function of Cr yields (see also [32] for a similar analysis) 

(Cr(^i)Cr(^2)Cr(4)) = i27rfd {^iki)b'j^i [P^Xki)PcM + perms] , (31) 

where P^^{k) = 2Tr'^VQ{k)/k^ and is a non-linear parameter given in this case by 



^^^"r (l + A-i)2 ' ^^^^ 



as follows from Eq. (26). Non-Gaussianities are thus significant when the curvaton decays well before 
it dominates, r <C 1. 

When A ^ 1 and the perturbations from inflation are negligible, one recovers the standard 
curvaton result [33] and is proportional to the much used local non-linear parameter /jvl 
defined by Cr = Cg + (3/5)/ArLCG' ~ (6/5)/iVL- However, in general is different from 

/jVL- Indeed, for other values of A, although only the curvaton contributes to the 3-point function, 
the 2-point function depends also on the initial inflaton fluctuation, which is a Gaussian random 
field, independent of the curvaton fluctuation. This differs from the original definition of /nl where 
only one Gaussian random field is present [34]. 

It is instructive to see how (32) depends on the curvaton expectation value during inflation, x*- 
Substituting the relation r ^ {x*l^pf' I \f^xh^ (valid in the limit r <C 1) [35], where is the 
decay rate of the curvaton, into the definition (30), one sees that A is proportional to x^, like r. One 
then finds that given in (32) reaches its maximal value 6^|(max) ~ e*/ ^/^^Jrrl^ for A ~ 1, 
i.e., for X* ~ i/r^/ {m^e^f)Mp. A significant non-Gaussianity is thus possible if 2> ■sjT^/m^. 
Note also that when r becomes small does not grow indefinitely as one would naively expect by 
considering /tvl — 5/(4r). Finally, in the limit r ^ 1 and A <C 1, where the inflaton contribution 
dominates the power spectrum, the expression (32) simplifies into 

3 £>f xi 



"^l^lT-iSrW' (33) 

After the analysis of the non-Gaussianities for the adiabatic perturbation, which essentially agrees 
with the discussion given in [32], let us now turn to entropy perturbations between CDM and 
radiation, 

<Sm = 3(Cm-Cr) , (34) 

which could be generated in the radiation era after the curvaton decay. If all the particle species 
are in full thermal equilibrium after the curvaton decays, with vanishing chemical potentials, then 
the primordial density perturbation must be adiabatic [6, 36] and we have Cm = Cr and hence 
5m = 0. However, if CDM remains decoupled from (part of) the radiation, the curvaton isocurvature 
perturbation may be converted into a residual isocurvature perturbation after the curvaton decays. 
We now consider two possibilities leading to a non-trivial isocurvature perturbation [6]. 

3.1 CDM created before curvaton decay 

If the CDM is created before the curvaton decay, then Cm = Cinf; which generates 



<5m = 3 (Cinf - Cr) = -tSg - - {- - 2r - r' ] S'q . (35) 
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This implies that the ratio between the isocurvature and adiabatic power spectra is given by 



Vr 1 + A 



_i • (36) 



This quantity is constrained to be small by the CMB data. In the case where the curvaton dominates 
the final C,, i.e. A S> 1, this scenario is thus ruled out and this case is often disregarded in the 
literature [7]. However, if the infiaton contribution is sufficiently important, A <C 1, such an entropy 
contribution is allowed. More specifically, the observational constraint on a, defined by Vs^jVc^^ = 
Q!/(l — a), is currently ao < 0.067 at 95% CL [37]. The subscript refers to the case where the 
entropy and adiabatic fluctuations are un- correlated^ which is appropriate here when A ^ 1. The 
non-Gaussianity of C is described by Eq. (33). Thus, it can become significant without violating the 
current bound on the presence of isocurvature component in the power spectrum. 

The amount of non-Gaussianity in the temperature fluctuations of the CMB anisotropies will 
depend both on ("r and S^. Thus, a complete study of these anisotropies would require the knowledge 
of the 3-point correlation properties of both variables. One can generalize Eq. (31) and define 

{X{kr)Y{k2)Z{k^)) = i2nf 5iJ:iki)b§r [PCrih)PcAk2) + perms] , (37) 

where X, Y, Z can be Cr or S^i- In this case all the 3-point functions have the same amplitude (up 
to numerical factors of —3), bff]^ = —Sbff^ = 9bf^^ = —27b'}^^, and equally contribute to the 
non-Gaussianity of the CMB temperature anisotropies. 



3.2 CDM created from curvaton decay 

The second possibility leading to non-trivial isocurvature perturbation is when the local matter 
density is produced solely from the local curvaton density (for instance, some fraction of the curvaton 
decays to produce CDM particles or the out-of-equilibrium curvaton decay generates the primordial 
baryon asymmetry) . Then we expect the matter density to be directly proportional to the curvaton 
density on the decay hypersurface 

Pcdme^^^-""^') = cp;,e3(^x-Cr) ^ (38) 

where c = (pm/Px) ^ hence to all orders = Cx- The matter isocurvature perturbation 

(34) is then given by 



("i -L fir- -L 2r'^ \ 
Sg - S'cj 

This implies that the ratio between the isocurvature and adiabatic power spectra is 

-Ps^ _ 9(1 -r)^ 

r2(l + A-i) ■ 



(39) 



(40) 



This quantity can be small, as required by observations, in two limiting cases. Either r is very close 
to 1 or A is very small. In the pure curvaton model (A S> 1), r is constrained to be very close to 1, 

9(1 - rf ~ a_i < 0.0037 (95% CL), (41) 
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using the constraint given in [37] for the totally anti- correlated case,^ and the non-linearity param- 
eters b-^J^^ defined in Eq. (37) involving entropy perturbations are suppressed by factors of (1 — r) 
with respect to , itself of order unity. 

If the inflaton dominates the linear perturbations, i.e. A <C 1, then the ratio (40) can be small 
even if the curvaton decays long before it dominates. For non-Gaussianity, the small r limit appears 
interesting because the amplitudes of the 3-point correlation functions b-^J^^ are related by 

bNL -^-bNL -^bj^L-:^b)^L^ (42) 

which shows that the amplitude of bff^ can be much larger than that of in this small r limit. 
However this situation is viable only if A satisfies the constraint 

V<: 9A 

^ ~ ^ ~ ao < 0.067 (95% CL) (r < 1, A < 1) , (43) 

using the constraint given in [37] for the un-correlated case. In this case is given by its limit in 
Eq. (33) and the expression for the non-Gaussianity of the isocurvature component, 



4f = -?(3 + 6r + 2.2)-Ji-445 (44) 



reduces to the simpler form 



2\r^J - * V 



X* 



^Xrlf^-TT ZT ^-326^ — , (45) 



where one recognizes the square of the isocurvature/adiabatic ratio Vs^/V(^^, which must be very 
small because of the constraint (43). 

Thus, even if the non-Gaussianity of the isocurvature component is much bigger than that of 
the adiabatic component, the constraint on the amplitude of the linear isocurvature perturbations 
(43) also constrains the magnitude of the bispectrum of isocurvature perturbations to be small. 
This follows simply from the fact that in this case the second order part of the matter isocurvature 
perturbation, in Eq. (39) is of order 5^, and if the linear part is constrained, then so is the 
non-linear part. 



4 Double quadratic inflation 

In this section, we consider double quadratic inflation [11], i.e., an inflationary phase driven by two 
massive and minimally coupled scalar fields described by the Lagrangian 

>C = -\{d^<l>f - \{d^xf - V{<l>,x) , V{4>,x) = \ml^^ + \mlx^ . (46) 

The adiabatic curvature perturbations generated by this model have been computed in [11, 38, 39] 
and their non-Gaussianity in [40, 41, 42] (see also [43]). This model can also generate isocurvature 

^Reference [37] uses the same convention as ours for the sign of the adiabatic and isocurvature perturba- 
tions. However, the cross-correlation is defined with opposite sign to {QSm), so that in the pure curvaton 
case adiabatic and isocurvature perturbations are referred to as being anti- correlated. 
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perturbations. At linear order these have been computed in [24] and it was first noticed in [44] that 
they can be correlated with the adiabatic ones. Here we extend these results at second order and we 
compute the 3-point correlation properties of the isocurvature perturbations and their correlation 
with the adiabatic perturbations. 

We start by computing the adiabatic curvature perturbation ^ at second order, using the 6N- 
formalism. The expression of the number of e-folds in terms of the scalar fields can be obtained from 
the slow-roll equations of motion, which read 

SiJ^Ml = im^</,2 ^ im^x^ , 3H<j> + ml<P = 0, 3Hx + mlx = 0, (47) 

and imply dN/dt = H = —{(pcj) + xx)/(2^p)- For a given scale, the number of e-folds between 
horizon crossing t^, and some subsequent time t is given by the expression 

N = ^^{<t>l + xl-<t>'-x') , (48) 

where (j) = (p{t) and x = x(^)- Furthermore, from the last two slow-roll equations in (47) one can 
derive 

= -, (49) 

X* 

where R = m^/nitf, and without loss of generality we take i? > 1. 

To compute we choose as final hypersurface at time t a uniform density hypersurface defined 
by the condition 

+ ^^ = C. (50) 

Then, Eqs. (49) and (50) uniquely fix the relation between (^, x) and (^*,x*)- Indeed, by combining 
these two equations we can derive 

R' xl + cl>' = c, 4,l = c, (51) 

which can be used to find the derivatives of and x with respect to 0^, and x*- These relations can 
then be employed to compute the first and second derivatives of N with respect to 0* and x* by 
differentiating Eq. (48). The calculation is reported in Appendix A. By using Eq. (15) up to second 
order in (5<^* and Sx*-, one obtains an expression for 

If- 11 



(i?2 - 1) 



^X* f-t _^X* \ _ A _ ^\ _^ X*9,x* f^X* j^2' 



X* V 2x* J 0* V 2^* y 2^ V X 



(52) 



where g = g{(f), x) is defined as g = <^x^ j {R'^X^ + (p"^)- This relation holds until the end of slow-roll 
infiation. 

We will consider the following scenario. Inflation is initially driven by the heavy fleld, X; which 
slow-rolls down the potential. Later, the heavy field becomes subdominant and then starts oscillat- 
ing, while the light field, (f>, drives inflation. In the last stage of slow-roll inflation, when x <S Mp, the 
coefficients g/Mp and X*9,x*/-^p "^^^y small and the second line of Eq. (52) becomes negligible. 
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The curvature perturbations ( thus becomes effectively constant, and since its value is unaffected 
by the subsequent stages of inflation and reheating, one finds the expression 

Cr = (^^*6(f>* + X*Sx* + + \^xl^ (53) 

for the adiabatic perturbation during the radiation era. 

Let us now focus on the isocurvature perturbation which can be produced, after inflation, in this 
type of model. To determine it, it is convenient to use the relative comoving curvature perturbations 
7?.0 and Tl^. According to Eq. (13), they are given by 

n^ = 6N-[^H^, n^ = SN-rH^. (54) 

The light field, (f), remains in slow-roll all the time during inflation.^ In order to compute TZ^ 
at second order we can use Eq. (14) for ipA = <P-i and expand N'^'^\ i.e., the number of e-folds from 
an initial fiat hypersurface at to a final uniform field (j) hypersurface, up to second order in 
and 5x*- To compute N^'t'^ we substitute Eq. (49) in Eq. (48) and impose that the final value of (p 
is a constant, (/> = C^. Using this condition and differentiating A^^"^) with respect to the initial field 
values (see Appendix A), Eq. (14) then yields 

where the explicit expression for C^j. is given in (53). At the end of inflation (j) dominates and 
reheats the universe. When ^ becomes dominant its comoving and uniform energy density curvature 
perturbations arc the same on large scales, = TZcj). Furthermore, when x^ <C Mj,, the second term 
on the right hand side of Eq. (55) becomes negligible, and = C,^. 

The evolution of the perturbation of the heavy field x is more complicated. During the (p- 
dominated slow-roll phase, a calculation similar to the one for (f) yields, replacing by x and by 
in Eq. (55), 

„ I _ ^ ^ ^ _ ^ _ _ iM ^ A^^^ 

/^Xlslow-roll - + ^2 ^^2 xl 2 <t^l ^ h xj ' ^ ^ 

where we have used that ^ dominates the Universe and thus H"^ = m'^^'^ /{6Mp). This expression is 
valid when x is subdominant and in slow-roll but cannot be used when x oscillates. It is convenient 
to use Eq. (54) to rewrite Eq. (56) in terms of the field fluctuation 6x on a constant total energy 
density hypersurface {SN = Q, 

K.u,-,„„=<-ri/?=<+^rf , (57) 

where we have used the property that H, which depends only on the slow-rolling is (spatially) 
constant on a constant total energy density hypersurface, since the latter coincides with a constant 



■^Wc assume that there is no intermediate dust-like phase between the heavy field dominated inflation and 
the light field dominated infiation. 
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(j) hypersurface when % is subdominant. This yields the non- linear relation between the isocurvature 
perturbation during slow-roll and the local value of the heavy field, 



^2 

^ = T^g5^(''^xlalow-roll-C) ^ j-gg^ 

which expanded up to second order yields 

?,H^ {5x I5x^ 



As for the curvaton, when x oscillates we can describe it as a non-relativistic fluid and use 
Eq. (18). Expanding this equation up to second order in the field fluctuation 6x^ we obtain the 
curvature perturbation C,^ during the oscillations, 

Cxiosc = Cr + o -7 - ■ (60) 



X 2 x' 

Now we can use the constancy of 6x/x valid up to second order to match Eqs. (59) and (60) and 
express Cxiosc in terms of 7?.^|siow-roii- Using Eq. (56), we find that the value of Cx after inflation is 

We assume that the light field decays into radiation which dominates the Universe after inflation, 
and that the heavy fleld decays into CDM when it oscillates. Then 

= 3(Cx - Cr) , (62) 

and 

This equation generalizes at second order the results of [24]. 

At this point it is useful to express the final curvature and entropy perturbations in terms of the 
instantaneous adiabatic and entropy perturbations during inflation (more precisely when the scale 
of interest exits the Hubble radius). The decomposition of two scalar field perturbations in terms 
of (instantaneous) adiabatic and entropy perturbations has been introduced at linear order in [12] 
and generalized at non-linear order in [15]. The general definitions are recalled in Appendix B. Here 
we give the expressions for the adiabatic perturbation 6a and the entropy perturbation Ss for the 
particular case of double quadratic inflation. To simplify the notation, let us deflne 

C0 = cos0 = -4>/i , se = sine = -R^x/^ , (64) 

with ^ = (0^ -I- R^x^y/'^. The angle 9 is simply the angle between the instantaneous direction of the 
fleld trajectory and the (f)-eLxis. At first order in perturbations we have 

Sa^^^ = ceSc/) + SffSx , Ss^^^ = cqSx — seScf) , (65) 

while the second order expression are given by 

2? 

C 2^ 
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Note that at second order the definition of 5a contains first order perturbations of 5s^^^ and vice 
verse. Indeed, as explained in [15], the adiabatic and entropy field decomposition is local and second 
order fluctuations will be sensitive to first order fluctuations of the angle 9, which can be re-expressed 
in terms of the field fluctuations 5(t^^^ and 5s^^\ 

Using these definitions and evaluating Eq. (52) at i = t*, we can rewrite (^*, i.e. the curvature 
perturbation on uniform density hypersurfaces at Hubble crossing. 



2Ml 



{eel + R'^sel) ^Ja, -I {I- R-\l - R^fceW:) 5al 



+ {R^ - l)ce^S0^{c0l + R "^801)50 Js^ . (68) 



The second order expression for C* contains also the first order entropy field perturbation. Indeed, 
this is the case also for its general slow-roll form (107) given in Appendix B, derived in [15]. 

The entropy field perturbation sources ^ at first and second order (see Eq. (102) in Appendix 
B). Using Eqs. (53) and (68), the final value of C is thus 



l-R~\ (, R'cet - se t 

Oil ^2 ^*(^0*Se* OS* — 



Cr = C*+ Oil ^2 ^*(^0*^0* — ^ — - — r^^^* ■ (69) 



Furthermore, we can rewrite the expression for Sm, Eq. (63), in terms of the adiabatic and entropy 
field perturbations. This reads 



C0*So*^* \ 2cg*S0*^^: 

Neglecting slow-roll corrections, the field perturbations 5(7* and 5s* are random fields with 2-point 
functions given by 

tt2 

{6a, {k)5a, {k')) = {5s* {k)5s* {k')) = i27rf 5{k + k')^ , {5a* {k)ds* {k')) = . (71) 

At lowest order in slow-roll, these fields are Gaussian. However, their 3-point correlation functions 
are non vanishing and have been computed in Appendix B. From the expression of C*, Eq. (68), and 
from the general definition of ( given in Appendix B, Eq. (107), one can sec that there are second 
order corrections proportional to 5a* and 5a*5s* so that (* has not exactly the same correlation 
properties as 5a*. However, these contributions are generically (for any slow-roll model) of the same 
order in slow-roll as the 3-point function of 5a* so that at lowest order in slow-roll C* oc 5a* is a 
Gaussian random field. 

A more convenient parametrization often used in the literature [11] is to rewrite the background 
values of the scalar fields in polar coordinates, 

= 2MpA//Ve^r/Vcosa , x = 2Mp a/A^^T/V sin a . (72) 

In terms of the angle a and of Ng — N*, the number of e- folds from Hubble crossing to the end of 
inflation, the power spectrum of C* is, using the linear term in Eq. (68), 

V -(N AM (l + ^^tan^«*)^ ^ V (73) 

^C* l^^e '^*\i^tein^a*){l + RHan^a*)\2TrMpJ ' ^ ' 
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Furthermore, instead of using (5s*, we find it useful to rewrite Eqs. (69) and (70) in terms of 
defined as having the same power spectrum of C*, i-e., 



2Ml 



(74) 



in analogy with the linear term of Eq. (68) . At leading order in slow-roll 5* is a Gaussian random 
field uncorrelated with C,^,. Finally, using this parametrization we obtain 



Cr = C* + 



(1 - Br) tana* 
1 + tan2 a* 



tan^ a* 



2 tana*(l + i?^tan^ a* 



and 



l + i?Maii- rt^ 



tana* 

where 77^,^ is a slow-roll parameter, 

V<l>4> 



r?^0 2 + + IR ^ I ail- a + R'^ t mii ^ q„ 
~2 i?2tana*(l -Fi?4tan2 a*) * 



1 + tan a* 



2(A/-e-iV*)(l + i22tan2a*) ' 



(75) 



(76) 



(77) 



As discussed in [44], adiabatic and entropy perturbations can be correlated at linear order but 
the correlation can be neglected when i^^tana* ^ 1 or tana* S> 1. Indeed, in this case Eqs. (75) 
and (76) yields, at linear order, = and <Sin oc 5*. However, these equations show that adiabatic 
and entropy perturbations are always correlated at second order, even when they are uncorrelated 
at linear order. In this particular model of two quadratic potential that we could treat analytically, 
non-linear terms are small, being suppressed by slow-roll. This leads to small non-Gaussianities in 
the adiabatic perturbation (cf. Ref. [40]) and also in the entropy perturbation. However, we do not 
expect this to be a generic feature of all inflation models. In particular, the coefficients in front 
of the S'^ terms in Eqs. (75) and (76) may be much larger in other models, which can lead to a 
non-vanishing 3-point correlation between the adiabatic and entropy perturbations. 



5 Conclusions 

We have calculated the second-order primordial curvature and isocurvature perturbations from two 
models of inflation in the early universe. In the first example of a mixed curvaton-inflaton model 
we assume the curvaton is an isocurvature field completely decoupled from the inflaton field driving 
inflation. In the second, double quadratic inflation model the two massive fields driving inflation are 
gravitationally coupled during slow-roll. 

The field perturbations at Hubble-exit during slow-roll inflation are effectively independent Gaus- 
sian random fields; their cross-correlation and non-linearities are suppressed by slow-roll parameters. 
However the coupled evolution on large scales after Hubble-exit can lead to cross-correlations at lin- 
ear order [44, 12] and we have calculated the correlations that arise at second-order. This can lead 
to non-vanishing bispectra for the primordial curvature and isocurvature perturbations and their 
cross-correlations . 

In both cases we find that the non-linear primordial curvature and isocurvature perturbations 
(15) and (16) can be given in terms of the adiabatic and entropy field perturbations at horizon-exit 
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during inflation, 




(78) 



(79) 



where N describes the expansion to a surface of uniform radiation density in the radiation dominated 
era and AA^ describes the difference between the expansion to hypersurfaccs of uniform radiation 
density and uniform matter density. Note that AA^ vanishes for adiabatic perturbations, i.e., when 
6s^ = 0. 

In the mixed curvaton-infiaton model we identify the inflaton field with the adiabatic pertur- 
bations during inflation, 5a*, and the curvaton field with entropy field perturbations, 6s*. The 
bracketed terms in Eq. (78) can be neglected at leading order in a slow-roll approximation. It is 
well known that C can have a significant non-Gaussianity when r <C 1 in the curvaton scenario (i.e., 
when A, the ratio between curvaton and inflaton contributions to the curvature power spectrum, 
is large). However, in this case the primordial isocurvature perturbations are constrained to be 
very small [7, 8]. We have shown that it is possible for a residual isocurvature perturbation to 
have a bispectrum which is much larger than that of the adiabatic perturbation if A is small, i.e., 
if the inflaton perturbation dominates the primordial curvature perturbation at first order, and if 
the CDM is produced by the curvaton decay. However observational constraints on the power spec- 
trum of isocurvature perturbations also constrains the bispectra to be small in this case. The most 
interesting situation is the scenario where the CDM is created before the curvaton decay, which is 
viable if the inflaton contribution dominates the linear power spectrum. In this case, it is possible 
to obtain a strong non-Gaussianity if e* S> ^JV^^Jm^ and we have found that the non-Gaussianity 
of the adiabatic component and of the isocurvature component are of the same order of magnitude. 

In double quadratic inflation the two canonical fields, (/> and x, are coupled gravitationally and 
the adiabatic and entropy field perturbations, ba* and (5s*, are in general a linear combination of 
the two canonical fields. We have obtained explicit expressions at second order relating the initial 
and final adiabatic and isocurvature perturbations. In this simple case of two uncoupled fields with 
quadratic potentials we find that the non-linearities are small, but this need not be the case in other 
models. Indeed, as shown in Appendix B, in general two field slow-roll inflation we expect that the 
terms in square brackets in Eq. (78) can be neglected at leading order in slow-roll. However, the 
remaining non-linear terms due to initial entropy perturbation need not be suppressed. It would be 
interesting to investigate non-Gaussianity of isocurvature perturbations in more general models. 

The bispectra for primordial curvature and isocurvature perturbations and their cross-correlations 
are presented for a general two-field model in Appendix C, and given at leading order in a slow-roll ex- 
pansion. These show that the non-zero primordial bispectra (in both curvature and isocurvature per- 
turbations and their cross-correlations) arise due to entropy field perturbations at Hubble-exit. Our 
results for the primordial curvature perturbation are consistent with the non-linear (5iV-formalism 
[4], derived in the large scale limit where the separate universes approach [45] is used to evaluate 
the perturbed expansion using the homogeneous background solutions. In single-field slow-roll in- 
flation the perturbations are adiabatic on large scales and the bispectrum is suppressed by slow-roll 
parameters [30, 31]. 

In multiple-field inflation there is the possibility of additional observational features which are 
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absent in single-field models. We have shown that this could include the contribution of isocurva- 
ture field perturbations to the bispectra (3-point functions) as well as the power spectra (2-point 
functions). It is interesting to note that, in principle, the isocurvature perturbations might domi- 
nate the primordial bispectrum in, for example, the CMB temperature anisotropics while remaining 
sub-dominant in the power spectrum. The non-Gaussian primordial perturbations predicted from 
Gaussian field perturbations during infiation are of a specific local form, but should be distinguishable 
from the local non-Gaussianity of the primordial curvature described by conventional /jvl parame- 
ter. The bispectrum of the isocurvature perturbations can be characterized by a new non-linearity 
parameter, and the cross-correlated bispectra yield additional parameters. But in curvaton models, 
for example, they are all determined by the single model parameter, r, and thus could provide a 
strong test of the curvaton scenario. The best constraints on specific models of non-Gaussianity are 
based on matched filtering techniques [46] . It will thus be important to develop optimized constraints 
for the non-Gaussianity of primordial isocurvature perturbations to obtain the optimal constraints 
on a wider range of theoretical models. 
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A Derivatives of N, N^^^ and A^*^^^ in double inflation 



If (f) and X are the values of the fields on a final uniform total density hypersurface, R^x^ + 0^ = C*; 
using Eq. (49) we obtain 
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The total number of e-folds in double inflation is given by 



(80) 




(81) 



which can be differentiate with respect to and x* to yield 
d(t) d(t) E? dx 



dx 



1 



(82) 



50* (j)*(t) ' dx* X*^ ' ^0* 0*X 
where g = (fP'x^ /{R^x^ + 0^). By differentiating N in Eq. (80) we get 



dx* 



9 , 



x*x 




(83) 



N 



X* 



A.* _ 



(84) 



2M|, 
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and 



N 

N 



V (p* (pi 



2Ml L 

^(1 - i?2)i?2% = AT^ . 



1 (^2_-l)Ml, 



2M|. 



X* 



1 



2M|, L 



1 + - 1) 



X* X% 



(85) 
(86) 
(87) 



If the final hypersurface is a uniform field ^ hypersurface, (f) = C^^ the number of e- folds (80) 
reads, 



1 



AMI 



2R^ 



which can be differentiated to give 



(88) 



(4') 



N 



2M|, 



X* 



2Mj, 



l + R 



* 



2X 



(89) 
(90) 



and 



N 



(0) 



(0) ^ 



iV('^) 
,x«x* 



2M|, 



AT 



{<t>) 



l-i?2^(l + 2i?2 



2M|. 



1-^ 

2 
* 



(91) 
(92) 
(93) 



Similar expressions can be found if we consider a final uniform field x hypersurface. In this case 



AtW 

1 A* 



AT' 



(x) _ 



X* 
2Ml 



2M| 



1 + 



(94) 
(95) 



and 



^,xix* 2Ml 



l-i?-2^(l + 2i?-2) , 



AT' 



(x) _ j^{x) _ 



P*X* 



R 



X* 

-2 a2 



N 



(x) 



:X* 

1 



Ml x*4>* ' 



2M|, 



p 



(96) 
(97) 
(98) 
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B Adiabatic and entropy field decomposition 



In this section we review the adiabatic and entropy decomposition approach at hnear and second 
order, during slow-roll inflation. It is possible to make a rotation in field space to identify the 
instantaneous adiabatic and entropy field perturbations along and orthogonal to the field trajectory. 
We will use the results of [15] generalizing the work of [12] (see [14] for an equivalent approach). At 
linear order, the adiabatic and entropy field perturbations are defined, respectively, as 

(5(7(^) = cos 9S(f) + sin 06x , Ss'^^'^ = - sin eScp + cos 95x , (99) 

where tan 6* = x/4> 6 is the time-dependent angle of the instantaneous rotation. At second order, 
we define 

5a = + ^ , Ss = ds^'^ - ^-^ ( 6s + %] , (100) 
2(7 a \ 2 I 

where a^ = <j)'^ + = -2HM],. 

The adiabatic curvature perturbation on uniform density hypersurfaces is defined at second order 

as 

(101) 



a (7 



-Scr] +l(-)5a + e^6s 
a I 2\ a I (7 



where bo is evaluated on a uniform flat hypersurface. The evolution of is sourced by first and 
second order perturbations of bs. On super-Hubble scales it reads^ 



2e&ds - {Vss + 4^^)5s^ + ^Ss6s 



a 



(102) 



where 6 = — F^/iT", with = — sinOV^^+cos 6Vy^ and Vgs = ^,4><t> sin^ ^~'^^,(t>x ^ ^™ ^+Kxx ^'^^^ ^■ 
The entropy field perturbation 5s evolves independently on super-Hubble scales and its evolution 
equation reads, 



5s + m5s + iVss + 3^^)(5s = ~5s^ -i(e + - 5s5s - \ ]-Vsss - ^^Vss - 9?^ 1 5s^ 

(7 a \ a 2 ) \2' '''' ^ 1 



(103) 

where = cos OV^j, + sin 6V^ and Vsss = —Vcpcpri) sin^ + SV^^./.^ cos 6 sin^ 9 — 3^^^^ cos^ sin ^ -|- 

Given these evolution equations we expect that their solutions can be written as 

C = C, + T^'Us,+T^'''^5sl, (104) 

6s = T^lhs, + T^^^dsl, (105) 

where T^^^^ and T^"^^^ are the first and second order transfer functions, for and 5s, respectively, and 
C*, (Js* are their initial conditions at Hubble exit. Equation (75) is an example of the general solution 
(104), while Eq. (76) is an example of (105) rewritten in terms of the CDM entropy perturbation 



^This equation corresponds to Eq. (221) of [15]. Note however that in vl and v2 of the arXiv and in the 
published version on JCAP, the last term inside the bracket is missing in this equation. We thank Sebastien 
Renaux-Petel and Gianmassimo Tasinato for pointing out this omission. 
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Let e be the standard slow-roll parameter, and let us define the mass slow-roll parameters rjij = 
Vij/ {SH'^) and rjas = {Vxx ~ '^h) Osm.O + ?70^(cos^ 9 - siv? 9), rjss = r]<P4> sin^ 9 — 2?7^^ cos 9 sin 9 + 

ri^^cos'^9. By using 9 = —Hrjas, which follows from the time derivative of tanS = x/(p ~ V^/V^^j), 
and 5s = —HrissSs, one can rewrite the second order definitions of da and ds, Eq. (100), and the 
definition of ^, Eq. (101), in terms of the slow-roll parameters, 

5a = 5a^^^ pi r]ssSs'^ , 5s = 5s^^^ + ( %,5s + ^5a] 5a , (106) 

and 



C 



5a 



1 



2eMp 



€ — 



2eMp 
^ 5a^ — ri„sSa5s 



(107) 



Note that in these three definitions, the non-linear terms on the right hand side are slow-roll sup- 
pressed with respect to the linear terms. If evaluated at Hubble crossing, when the slow-roll pa- 
rameters are small, they lead to contributions to the intrinsic non-Gaussianities of (5(7*, 5s*, and (^* 
which are small. 

For completeness, we compute here the 3-point correlation functions of 5(7* and 5s*. Prom the 
results of [26], namely 



(5(^i(fci)5(^*^(fc2)5(^f (4)) = i2TrfSiEiki 



Ht 



16M| 



perms 



HMikf 



with 



M{ki,k2,k3) = -kikl-4^ + hf+ ,2 

kt 2 Kt 



^21,2 
"'2 "'3 



(^2 - k^), kt = ki + k2 + kz 



(108) 



(109) 



one can compute the 3-point correlators of 5cr(^) and 5s^^\ simply by using the change of basis in 
field space (99). We find that {5a'i^ 5a'^}ha''}^) is the same as for a single scalar field [30], whereas 
(5(T*^'*5(7*^^5sl^^) and (5sl^^5sl^'*5sl^^) vanish, simply because s = 0. We also find 



(5sl^^(^i)5s:;'^(A:2)5c7r'(fc3)) = (27r)^5(Eife 



,(1) 



(1)/ 



T.Ht (kl-kz{kl + kl) 



2l2 



Uikf 



ktllikf 



(110) 



8V2Mp 

and similar expressions for {5s^^\ki)Sai^\k2)5s':^\k3)) and {5ai^\ki)5s^^\k2)5s^*Hk3)) by relabel- 
ing appropriately the ki appearing on the right hand side of (110). 

Taking also into account the second order parts of the adiabatic and entropy perturbations given 
in Eq. (106), we eventually find for the full 3-point functions 

-• , ^ 1 „ , -> . ^ffZHt I ^'i.kf Yjj^ikik^ Yji-^jki k^ \ 

{27rf 5{^iki 



(5(7* (A;i)5(7* (A;2)5c7* {ks)) 

(5(7* (fcl)5(7* (fc2)5s* (^3)) 
(5s*(^i)5s*(^2)5(7*(^3)) 

(5s*(^i)5s*(fe2)5s*(^3)) 



i2Trf 5i^iki 



8V2Mp 

. ^crs*H^, 



Ilikl 



ktUikf 



1 



8V2^Mp kfk^ 



{2nf5{T.ik 



Ht 



. 



8V2e;Mp 



3\ n 



^3 ^3(^1 ^2 



> ^l'^2 

'hliik^i 



(111) 

(112) 

(113) 
(114) 



This shows that the intrinsic non-Gaussianities of 5(7*, 5s* and, as a consequence of Eq. (107), of C* 
are all small for slow-roll models. 
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C Primordial power spectra and bispectra 

At first-order the expressions in Eqs. (78) and (79) for the primordial curvature and isocurvature 
perturbations in terms of the field perturbations during inflation give the power spectra of the 
primordial perturbations at leading order 

(Cr(^l)Cr(^2)) = {27r)H^{h + {N^aPAh) +2N sCsaih) + Nf.Psih)} , (115) 
^(Cr(^l)<5m(fc2)) = {2TT)^S^{kl + ^2) {NsANsPsih) + N^ANsCsaih)} , (116) 

l{Sra{ki)Sm{k2)) = i2TTf6^ih + k2)ANf^Ps{kl), (117) 

where at horizon-crossing during inflation we have 

{6a,{ki)6a,{k2)) = i27rf6{ki + k2)Pa{ki), (118) 
{5s,{ki)6s4k2)) = i27rf5{ki + k2)Psiki), (119) 
{Sa,ih)Ss,{k2)) = {27rf6{ki + k2)Csa{ki). (120) 

At second-order Eqs. (78) and (79) give the leading order bispectra for the primordial curvature 
and isocurvature perturbations and their correlations, 

(Cr(fcl)Cr(^2)Cr(4)) = NjNjNK{S'pi{ki)Svi{k2)S<pfih)) 

+ {2TrfS{^iki)NiNjNKL[C^^iki)C^^{k2) + 2 perms] , (121) 

l{UhKrik2)S^ih)) = NiNjANK{S<piiki)S<piik2)Scpfih)) 



+ 
+ 



(27r)3(5(E,^,) { NiNjANKLC'^'iki ) C"^^ (^2 ) 
NuNKANLlC-'^'iki) + C"'^(A;2)]C^^(A;3)} , (122) 



l{Sra{ki)Sra{k2)Uh)) = ANjANjNK{Sdiki)6ipi(k2)Sip^{h)) 
y 



+ 
+ 



{27rf5iJ:iki)[ANiANjNKLC'''iki)CJ'^ik2) 
ANjjANKNL[C'''{k,) + C-^-^(A;2)]C^^(fc3)} , (123) 

^{Sm{ki)S^{k2)S^{h)) = ANjANjANK{S^i{ki)5^i{k2)5^Hh)) 

+ {2Trf6{J:iki)ANiANjANKL[C^^{ki)C^^{k2) + 2 perms](,124) 

where C^^k) = Pi{k). 

At leading order in a slow-roll expansion the adiabatic and entropy field perturbations are inde- 
pendent Gaussian random fields, with [47, 48] 

Ps{k) ~ P„{k) ~ P,(fe) = ^ , CsAk) ~ , (125) 
and the primordial bispectra simplify considerably in the slow-roll limit, where we drop the terms 
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in brackets in Eqs. (78) and (79), to give 



iCrikMkXrih)) 

^(Cr(^l)Cr(^2)<Sm(^3)) 



{2TrfS{^iki){NsNssANs [P*{ki) + P*ik2)] P*{k3) 
+N^,ANssP*iki)P.ik2)] , 



{27Tfd{J:ik^NlNss [P*{ki)P^{k2) + 2 perms] , 



(126) 



(127) 



\{Sr.{k)Su.{k2)Uk3)) 



(27r)3(5(Eifci){iV,,AiV,,AiV,,, [P,{ki) + P,{k2)] P.{k^) 
+NssAN^sP*{ki)P*ik2)}, 



(128) 



■^{Sm{kl)Smik2)Sra{k3)) 



i2TTf6i^iki)ANf,AN,ss [P.{ki)P,{k2) + 2 perms] . 



(129) 
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